Spin-fluctuation mechanism of superconductivity in cuprates 
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The theory of superconductivity within the t-J model, as relevant for cuprates, is developed. It is based 
on the equations of motion for projected fermionic operators and the mode-coupling approximation for the 
self-energy matrix. The dynamical spin susceptibility at various doping is considered as an input, extracted 
from experiments. The analysis shows that the superconductivity onset is dominated by the spin-fluctuation 
contribution. We show that T c is limited by the spin-fluctuation scale T and shows a pronounced dependence on 
the next-nearest-neighbor hopping t'. The latter can offer an explanation for the variation of T c among different 
families of cuprates. 

PACS numbers: 71.27.-Ha, 74.72.-h, 74.20.Mn 
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Since the discovery of high-temperature superconductivity 
(SC) in cuprates the mechanism of SC in these compounds 
represents one of the central open questions in the solid state 
theory. The role of strong correlations and the antiferromag- 
netic (AFM) state of the reference insulating undoped com- 
pound has been recognized very early 1 1]. Still, up to date 
there is no general consensus whether ingredients as embod- 
ied within the prototype single-band models of strongly corre- 
lated electrons are sufficient to explain the onset of high T c , or 
in addition other degrees of freedom, as e.g. phonons, should 
be invoked. As the basis of our study, we assume the simplest 
t-J model [2], allowing besides the nearest-neighbor (NN) 
hopping t also for the next-nearest-neighbor (NNN) hopping 
t' term. The latter model, as well as the Hubbard model |3], 
in the strong correlation limit U ^> t closely related to the t-J 
model, have been considered by numerous authors to address 
the existence of SC due to strong correlations alone. Within 
the parent resonating-valence-bond (RVB) theory |[ll 0] and 
slave-boson approaches to the t-J model |4] the SC emerges 
due to the condensation of singlet pairs, induced by the ex- 
change interaction J. An alternative view on strong correla- 
tions has been that AFM spin fluctuations, becoming particu- 
larly longer-ranged and soft at low hole doping, represent the 
relevant low-energy bosonic excitations mediating the attrac- 
tive interaction between quasiparticles (QP) and induce the 
d-wave SC pairing. The latter scenario has been mainly fol- 
lowed in the planar Hubbard model 1 3] and in the phenomeno- 
logical spin-fermion model @|. Recent numerical studies of 
the planar t-J model using the variational quantum Monte 
Carlo approach as well as of the Hubbard model using 
cluster dynamical mean-field approximation |7], seem to con- 
firm the stability of the <i-wave SC as the ground state at inter- 
mediate hole doping. 

The t-J model is nonperturbative by construction, so it is 
hard to design for it a trustworthy analytical method. One 
approach is to use the method of equations of motion (EQM) 
to derive an effective coupling between fermionic QP and spin 
fluctuations 1 8 ] . The latter method has been employed to eval- 
uate the self energy and anomalous properties of the spectral 
function 1 8 . l9t 1 1 Oil . in particular the appearance of the pseudo- 



gap and the effective truncation of the Fermi surface (FS) at 
low hole doping fioll . The analysis has been extended to the 
study of the SC pairing while an analogous approach has 
been also applied to the Hubbard model II ill . 

In the following we adopt the formalism of the EQM and 
the resulting Eliashberg equations within the simplest mode- 
coupling approximation |*J 1 1(11 . Equations involve the dy- 
namical spin susceptibility which we consider as an input 
taken from the inelastic-neutron-scattering (INS) and NMR- 
relaxation experiments in cuprates. The analysis of these ex- 
periments 111 211 reveals that in the metallic state the AFM stag- 
gered susceptibility is strongly enhanced at the crossover from 
the overdoped (OD) regime to optimum (OP) doping and is 
increasing further in underdoped (UD) cuprates, while at the 
same time the corresponding spin-fluctuation energy scale is 
becoming very soft. Direct evidence for the latter is the ap- 
pearance of the resonant magnetic mode within the 
SC phase indicating that the AFM paramagnon mode can be- 
come even lower than the SC gap. These facts give a renewed 
support to the scenario that the spin fluctuations in cuprates 
represent the lowest bosonic mode relevant for the mechanism 
of the d-wave SC pairing. One of the novel results of our ap- 
proach is the importance of the NNN hopping t', consistent 
with the emerging evidence for the variation of T c among dif- 
ferent families of cuprates 1 15]. 

We consider the t-J model 



H 



<ij) 



(1) 



including both the NN hopping iy = t and the NNN hopping 



tij = t . The projection in fermionic operators, Ci S = (1 — 
n,i t -s)cis leads to a nontrivial EQM, which can be in the k 
basis represented as 



[c ks ,if] = [(l + C/l )| 



J(l - c h )]c ks + 



(2) 



1 



where Ch is the hole concentration and mtq is the effective 
spin-fermion coupling = 2J7 q + e k _ q , while e£ = 



2 



— 4^7k — 4t'7£. is the bare band dispersion on a square lattice. 
As argued in Ref. 1 10] we use the symmetrized coupling in or- 
der to keep a similarity with the spin-fermion phenomenology 



TOkq = 2J7 q + i(e£ 



k-q 



el). 



(3) 



EQM, Eq. (0, are used to derive the approximation for the 
Green's function (GF) matrix Gks(^) = ((^ksl for 
the spinor "fks = (ck, S ; e_ k _ s ). We follow the method, as 
applied to the normal state (NS) GF by present authors jjlfioll. 
and generalized to the SC pairing by Plakida and Oudovenko 
|0]. In general, we can represent the GF matrix in the form 

GksM" 1 = ~[WT - Cks + l iT i - S ksM], (4) 

a 

where a = J2i({^is,c\ s } + )/N = (1 + c^)/2 is the nor- 
malization factor, /i is the chemical potential and £ks the fre- 
quency matrix, 



6» = x<{[*k», 



(5) 



Eq. generates a renormalized band £k = = £ — 
4f7i*7k — 477 2 i'7 k and the mean-field (MF) SC gap 

A k = C k ? = --^5>k- q <c_ q ,- s c q , s >. (6) 
q 

To evaluate Ek s (w) we use the lowest-order mode-coupling 
approximation, analogous to the treatment of the SC in the 
spin-fermion model [5], introduced in the t-J model for the 
NS GF [ni and extended to the analysis of the SC state 
Taking into account EQM, Eq. (0, and by decoupling 
fermionic and bosonic degrees of freedom, one gets 



v ll(12),. s 
S k S 



-1 

Na(3 



££<«> m ) 



q m 



[3x q (iw„ - iui m ) ± -X°Jiu n - iuo n 



(7) 



where iu n — in(2n + l)/[3 and Xq(w),x q ( w ) are dynami- 
cal spin and charge susceptibilities, respectively. In cuprates 
charge fluctuations seem to be much less pronounced, so we 
furtheron keep only the spin contribution. 

In order to analyze the low-energy behavior in the NS and 
in the SC state, we use the QP approximation for the spectral 
function matrix 



2£ k 



( wt - e k T 3 - A ks Ti ) [5 ( w - i? k ) - <S + i?k ) ] , 

(8) 

where £k = (e k + A^) 1 / 2 , while NS parameters, i.e., the 
QP weight Zk and the QP energy ek, are determined from 
Gks(^ ~ 0), Eq. @. The renormalized SC gap is 



v 12 
^ks 



(0)]. 



(9) 



It follows from Eq. (@J that G k 2 (ia;„) ~ -a2 , k A ks /(w 2 + 
E%). By defining the normalized frequency dependence 
Fq(iu)i) = Xq(* w 0/Xq' an( l rewriting the MF gap, Eq. ©, 
in terms of the spectral function, Eq. l|8}, we can display the 
gap equation in a more familiar form, 



Aks = ^£[ 4J ^k-q-3TOk,k-qXk-qCq,k-q] 



q 

(Z0^A qs /2£ q )th(/3£: q /2) 



(10) 



where Ck q = ik q (^n ~ 0)// k plays the role of the cutoff 
function with 



Ikq{iv n ) = 4 £ Fq(iuj n - ioj m ) 2 ^ 



(11) 



and /£ = tb.(j3Ek/2)/(2Ek). Eq. <[To} represents the BCS- 
like expression which we use furtheron to evaluate T c , as well 
to discuss the SC gap A q (T = 0). To proceed we need the in- 
put of two kinds: a) the dynamical spin susceptibility % q (u;), 
and b) the NS QP properties Z^, e^. 

The INS experiments show that within the NS the low-w 
spin dynamics at q ~ Q is generally overdamped in the whole 
doping (but paramagnetic) regime [14]. Hence we assume 
X q (w) of the form 



BqUJ 



r 2 : 

q 



Fq(iui) 



(12) 



Following the recent memory-function analysis 1 16] B q = 
X q T q should be quite independent of q = q — Q. We choose 



the variation as T r 



r Q (i 



wq 2 /k 2 ) 2 consistent with the 



INS observation of faster than Lorentzian fall-off of x q (w) vs. 
q 1 14]. w ~ 0.42 in order that k represents the usual inverse 
AFM correlation length. 

Consequently, we end up with parameters XqjFq,k, 
which are dependent on Ch, but in general as well vary with 
T. Although one can attempt to calculate them using the anal- 
ogous framework llal . we use here the experimental input for 
cuprates. We refer to results of the recent analysis 11211 . where 
NMR T2G relaxation and INS data were used to extract k, 
Xq (T) and Tq (T) for various cuprates, ranging from the UD 
to the OD regime. For comparison with the t-J model, we use 
usual parameters t = 400 meV ~ 4500 K, J = 0.3i. At least 
for UD cuprates, quite consistent estimates for Xq,Tq can 
be obtained also directly from the INS spectra 11411 . For UD, 
OP and OD regime, i.e., c h = 0.12, 0.17, 0.22, respectively, 
we use furtheron the following values: Xq^ = 15.0, 4.0, 1.0, 
T° Q /t = 0.03,0.1,0.18 (appropriate at low T), and k = 
0.5, 1.0, 1.2. It is evident, that in the UD regime the energy 
scale Tq becomes very small (and consequently Xq oc 1/Tq 
large, in spite of modest k 111 ), supported by a pronounced 
resonance mode lfl4ll . We take into account also the T depen- 
dence, i.e., Tq(T) ~ Tq + T 1 12], being significant only in 
the UD regime. 

For the NS Ak(ui) and corresponding Z^.e^ we solve 
Eq. for S^ 1 = S k as in Ref. [1(3], with the input for XqM 
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as described above. Since our present aim is on the mech- 
anism of the SC, we do not perform the full self-consistent 
calculation of Ekfw), but rather simplify it as done in the pre- 
vious study lHoll . As evident from EQM, Eq. (0, the effective 
spin-fermion coupling mtq is strong (of the order of t) gener- 
ating a strong incoherent component of the spectral function, 
leading to an overall decrease of the QP weight Z < 1 and 
the QP dispersion with renormalized 771 , 772 < 1. Soft AFM 
fluctuations with q ~ Q lead through Eq. @ to an additional 
reduction of Zk, which is k-dependent. A pseudogap appears 
along the AFM zone boundary and the FS is effectively trun- 
cated in the UD regime with Zk F -C 1 near the saddle points 
(tt, 0) (in the antinodal part of the FS) 1101 . We fix fx with the 
FS volume corresponding to band filling 1 — c/j. 

We first comment general properties of the gap equation, 
Eq. JlOi . Close to half-filling and for x q peaked at q ~ Q 
both terms favor the d x i_ y i SC. The MF-part A k , Eq. (jfji, 
involves only J which induces a nonretarded local attraction, 
playing the major role in the RVB theories In contrast, 

the spin-fluctuation part represents a retarded interaction due 
to the cutoff function Ck q determined by Tk- q . The largest 
contribution to the SC pairing naturally arises from the antin- 
odal part of the FS. Meanwhile, in the same region of the FS 
also Zk is smallest, reducing the pairing strength in particular 
in the UD regime. Our analysis is also based on the lowest 
order mode-coupling treatment of the SC pairing as well as 
of the QP properties near the FS. Taking this into account, 
one can question the relative role of the hopping parameters 
t, t' and the exchange J in the coupling, Eq. While our 
derivation within the t-J model is straightforward, an analo- 
gous analysis within the Hubbard model using the projections 
to the lower and the upper Hubbard band, respectively, would 
not yield the J term within the lowest order since J oc t 2 . This 
stimulates us to investigate in the following also separately the 
role of J term in Eq. i ll Oi l, both through the MF term, Eq. (|6jl, 
and the coupling mkq, Eq. 0. 

Let us turn to results for the NS spectral properties and con- 
sequently T c . As discussed above we assume that the effec- 
tive band is renormalized due to the incoherent hopping 1 10], 
with: 771 = 772 = 0.5, Z — 0.7, whereas the coupling to 
low-energy AFM fluctuations, Eq. 0, leads to an additional 
QP renormalization. For fixed t'/t = —0.3, we present in 
Fig. 1 results for the variation of the Z^ in the Brillouin zone 
for two sets of parameters, representing the UD and the OD 
regime, respectively. The location of the renormalized FS is 
also presented in Fig. 1. While the coupling to AFM fluctu- 
ations partly changes the shape of the FS, more pronounced 
effect is on the QP weight. It is evident from Fig. 1 that Z^ is 
reduced along the AFM zone boundary away from the nodal 
points. Particularly strong renormalization Z^ <C 1 happens 
in the UD case, leading to an effective truncation of the FS 
away from nodal points fioll . 

NS results for Z^, are used as an input for the solution 
of the gap equation, Eq. dlOi . as presented in Fig. 2. For the 
same t'/t = -0.3 we calculate T c /t for c h = 0.12,0.17,0.22. 
Besides the result a) of Eq. dlOt (full line in Fig. 2) we present 
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Figure 1: QP weight evaluated for t'/t = —0.3 for parameters 
corresponding to Ch = 0.12 and cj, = 0.22, respectively. White line 
represents the location of the FS. 



also two alternatives: b) the solution of Eq. dlOt without the 
MF term, and c) the result with TOk q without the J term and 
omitted MF term. In the latter case, we used as input NS QP 
parameters, recalculated with correspondingly modified rrikq- 
From Fig. 2 it is evident that the spin-fluctuation contri- 
bution is dominant over the MF term. When discussing the 
role of the J term in the coupling, Eq. Q, we note that in 
the most relevant region, i.e., along the AFM zone boundary 
TTikQ = 2 J — At' cos 2 k x . Thus, for hole doped cuprates, 
t' < and J terms enhance each other in the coupling, and 
neglecting J in ?7ikq reduces T c , although at the same time 
relevant Zk is enhanced. 

0.04 1—1 ■ 1 . 1 i 1 




1— ' ■ 1 1 ' 1 1 

0.1 0.15 0.2 0.25 

c k 

Figure 2: T c /t vs. doping for t'/t = —0.3, calculated for various 
versions of Eq. 1 1 ( )l : a) full result (full line), b) with neglected MF 
term (dashed line), and c) in addition to b) modified 77ik q without the 
J term (dotted line). 

Finally, in Fig. 3 we present results, as obtained for fixed 
OP c h = 0.17, but different t'/t < 0, as relevant for hole- 
doped cuprates 1 15]. As expected the dependence on t' is pro- 
nounced, since the latter enters directly the coupling ?rikQ. It 
is instructive to find an approximate BCS-like formula which 
simulates our results. The latter involves the characteristic 
cut-off energy Tq, while other relevant quantities are the elec- 
tron density of states A/q and Z m being the minimum Z^ on 
the FS (in the antinodal point). Then, we get a reasonable fit 
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0.1 0.2 0.3 0.4 

-fit 

Figure 3: T c /t vs. —t'/t for fixed 'optimum' doping c;, = 0.17 and 
different versions ofEq. (TO), as in Fie. 2. 

to our numerical results with the expression, 

T c ~ 0.5r Q c -2/(AroV e// ) ; (13) 

where the effective interaction is given by V e f / = 3Z m (2J — 
4£') 2 xq. Our numerical analysis suggests that the main t'- 
dependence of T c originates in the coupling rh^q, not in 
AfoZ m , while the main ct -dependence comes from xq and 
Tq. Then, Eq. (1131 implies that optimum doping, where T c 
is maximum, increases with —t'/t. For parameters used in 
Fig. 1, e.g., c opt = 0.13 + 0.12(-i'/i). 

In conclusion, let us comment on the relevance of the 
present method and results to cuprates. Our starting point is 
the model, Eq. Q, where strong correlations are explicitly 
taken into account via the projected fermionic operators. In 
this respect the derivation crucially differs from the analysis 
of the phenomenological spin-fermion model |5]. Neverthe- 
less, in the latter approach the resulting gap equation, Eq. i ll Oi l, 
looks similar but involves a constant effective coupling. In 
contrast, our m^, Eq. (0, is evidently k, q- dependent. In 
particular, in the most relevant region, i.e., along the AFM 
zone boundary, mkQ depends only on t' and J, but not on t. 
This explains our novel result, i.e., a pronounced dependence 
of T c on t' which is also consistent with the evidence from 
different families of cuprates 1 15]. One can give a plausible 
explanation of this effect. In contrast to NN hopping t, the 
NNN t' represents the hopping within the same AFM sublat- 
tice, consequently in a double unit cell fermions couple di- 
rectly to low-frequency AFM paramagnons, analogous to the 
case of FM fluctuations generating superfluidity in 3 He [vfa . 

It is evident from our analysis, that actual values of T c 
are quite sensitive on input parameters and NS properties. 
Since we employ the lowest-order mode-coupling approxi- 
mation in a regime without a small parameter, one can ex- 
pect only a qualitatively correct behavior. Still, calculated T c 
are in a reasonable range of values in cuprates. We also note 
that rather modest 'optimum' T c value within presented spin- 
fluctuation scenario emerge due to two competing effects in 
Eqs. < l 1 01 . ( 1 1 3I > : large mtq and xq enhance pairing, while at 
the same time through a reduced and cutoff Tq they limit 



We do not discuss in detail results for Ak in the SC 
phase, where we obtain the expected d x 2_ y 2 form with At ~ 
A (cos 2 fc ;E - cos 2 fc y )/2, with A (T = 0) ~ r)T c and 
i] ~ 2.5. However, we observe that in the UD regime the 
effective coherence length £ ~ vf/Aq(T = 0) becomes very 
short. I.e., with vp taken as the average velocity over the re- 
gion k at the antinodal part of the FS we get £ ranging from 
£ = 4.4 in the OD case, to £ = 1.3 in the UD example. In 
the latter case, SC pairs are quite local and the BCS-like ap- 
proximation without phase fluctuations, Eqs. J 1 Ofc . J 1 3i . over- 
estimates T c . Starting from this side, a more local approach 
would be desirable. 

It should also be noted that in the UD regime we are deal- 
ing with the strong coupling SC. Namely, we observe that 
AfoV e ff shows a pronounced increase at low doping mainly 
due to large xq. Then it follows from Eq. Jl 3I > that T c is 
limited and determined by Tq. At the same time, INS exper- 
iments 1 14] reveal that in th UD cuprates the resonant peak at 
uj ~ u> r takes the dominant part of intensity of q ~ Q mode 
which becomes underdamped possibly even for T > T c . Thus 
it is tempting to relate Tq to u> r (for more extensive discus- 
sion see Ref. flfilo and in the UD regime to claim T c ~ au> r , 
indeed observed in cuprates fl4Tl with a ~ 0.26. However, 
additional work is needed to accommodate properly an under- 
damped mode in our analysis. 
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